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1. Introduction 


The two-pion production in nucleon-nucleon (NN) collisions 
is a rich source of information about the baryon excitation spec¬ 
trum and the baryon-baryon interactions. In addition to the ex¬ 
citation of a resonance decaying into two pions, which can also 
be studied in the nN nnN HI] and jN — > imN 12] reac¬ 
tions, the simultaneous excitation of two baryons can be inves¬ 
tigated in the NN reactions. By giving access to single and 
double baryon excitation processes, which both play an impor¬ 
tant role in the NN dynamics in the few GeV energy range and 
contribute significantly to meson and dilepton production, the 
two-pion production appears as a key process towards a bet¬ 
ter understanding of hadronic processes. In comparison to the 
one-pion decay mode, it presents a different selectivity with re¬ 
spect to the various resonances. In particular, the excitation 
of baryonic resonances coupled to the p meson can be studied 
with the two pions in the isospin 1 channel. This is of utmost 
interest for a better understanding of the dilepton production 
in nucleon-nucleon reactions, where these couplings manifest 
clearly llii , and also in nucleon matter due to the expected 
modifications of the p meson spectral functions @]. Finally, 
the comparison of two-pion production in pp and np channels 
could shed some light on the origin of the surprisingly large 
isospin dependence of the dilepton emission observed by the 
HADES experiment 12. In particular, the p production mecha¬ 
nism via AA final state interaction, which does not contribute in 
the pp channel, was recently proposed as an explanation for the 
different dilepton yield measured in pp and pn channels |2. It 
is therefore important to check the description of the double A 
process in the two-pion production channels. 

Additionally, following the intriguing results obtained by the 
WAS A collaboration in the double pionic fusion reactions, a re¬ 
newed interest on the study of the two-pion production in NN 
collisions was sparked, in order to check the possible contribu¬ 
tion of a dibaryon resonance M.M- 

The answer to all these open questions requires system¬ 
atic two-pion production measurements both in proton-proton 
and neutron-proton collisions. Concerning proton-proton col¬ 
lisions, a significant amount of data has been accumulated for 


various two-pion final channels in bubble chamber experiments 

for proton 

incident energies from the threshold up to 2.85 GeV. Precise 
differential cross-sections have also been obtained recently at 
CELSIUS and COSY up to 1.4 GeV 

3lLl32Ll33Ll34l] . with an emphasis on the production. The 


data base for the pn reaction from the bubble chamber experi¬ 
ments is even more scarce j 3 18, 3 21, 35). Very recently, 
however, precise measurements of total and differential cross 
sections for the np —> ppir^n^ and np —> npuPiiP became avail¬ 
able from WAS A at COSY at neutron energies from 1.075 to 
1.36 GeV |3. 22I- In the np —> npn^n" channel, differen¬ 
tial cross-sections are also known from Dubna measurements 
138. 3] , covering the beam incident energy range from 0.624 
to 4.346 GeV. 


10125 Torino, Italy. 


Since the chiral perturbation theory calculations for two-pion 
prod uction in NN collisions are available only near threshold 
||40|] , several phenomenologic al models have been suggested for 
the analysis of the double pion production in NN collisions in 
the GeV energy range. The first theoretical developments re¬ 
lated to the two-pion production were based on the one-pion 
exchange (OPE) model The reggeized n exchange model 
(OPER) ll42[ 1431] . which uses the partial wave analysis results 
for 7tN elastic scattering ll42 . constitutes its most recent and 
most elaborate modification. Lagrangian models were also in¬ 
troduced. The Valencia model by Alvarez-Ruso et al. ll42 was 
first developped. It aimed at a description of NN collisions at 
energies lower than 1.4GeV and included A( 1440) and A(1232) 
excitations. The Cao et al. model developped after the 
publication of new data at COSY and CELSIUS, has a larger 
range of applicability due to the inclusion of resonances with 
mass up to 1.72 GeV. Both models qualitatively reproduce the 
very fast increase of the cross section above threshold in the dif¬ 
ferent two-pion production channels and predict the dominance 
of two processes above 1 GeV: the excitation of the N(1440) 
resonance and subsequent decay into An or Ncr and the double 
A excitation. 

However, both models ] 3 have failed to reproduce the 
n'^n^ spectra for the ^ ppn^n^ reaction at beam ener¬ 
gies above 1.0 GeV 132 


3411 . which motivated the develop¬ 


ment of the so-called ’’modified Valencia model” 0321] . provid¬ 
ing a much improved description of these data. This new model 
has been used by the WASA collaboration for the interpreta¬ 
tion of the double pionic fusion reactions, after some additional 
changes to take into account the deuteron formation. How¬ 
ever, the observed resonant behavior of the cross section of the 
pn dn^n^ ]S[I2. associated with a structure at low n^n^ in¬ 
variant mass (the so-called ABC effect) could not be explained 
by such an approach and were interpreted as being due to a 
dibaryon resonance in the 1=0 NN channel, with a mass of 2.37 
GeV/c^ and a width of 70 MeV. This hypothesis was further 
supported by the isospin decomposition of the pn dnn reac¬ 
tion ii. The latter provided a consistent description of both 
1=0 and 1 = 1 channels by taking into account the resonant con¬ 
tribution in addition to the conventional t-channel processes de¬ 
scribed by the ’’modified Valencia model”. Even more recently, 
the pn ppn^n^ 13] reaction was also consistently described 
with the same model. The accuracy of d* resonance hypothesis 
is also supported by the SAID partial wave analysis based on 
new polarized np scattering data. In this analysis a resonance 
pole in the G 3 coupled pmtial waves at (2380+10-140+5) 
MeV have been discovered ]^, 3l . 

Considering the impact that the discovery of a dibaryon res¬ 
onance could have, this systematic study must be pursued. It 
is indeed important to provide, possibly in independent exper¬ 
iments, constraints for all possible channels of the pn reac¬ 
tion, where the resonance is expected to contribute, in partic¬ 
ular n^n^ or n^n^ production with an unbound pn pair, but also 
for pp channels in order to check unambiguously the consistent 
description of the conventional processes. 

The experiments with HADES 13 are particularly well 
suited to answer such a request, since both pp and pn exper- 
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iments were studied in the relevant energy range. The good 
capacities of HADES for hadron identification were already 
shown in the study of one-pion and one-eta (measured by its 
three pion decay) production channels iBEl of the pp reac¬ 
tion at different energies. Beyond the search for the dibaryon 
resonance, the study of the two-pion production process by the 
HADES collaboration is also motivated by the connection with 
the dilepton production, as mentioned above. 

As a hrst step in this program, we report here on the analysis 
of the tagged quasi-free pn —> pnji^ji" reaction from experi¬ 
ments using a deuteron beam of 1.25 GeV/nucleon and present 
precise total and differential cross-sections. According to previ¬ 
ous estimates 114, 1_8, 21 ], this channel is dominated by the 

isospin-0 contributions, in which the resonance should reveal. 
In the present analysis, the results were averaged over the avail¬ 
able range of np center-of-mass energies. Our strategy consists 
in comparing the various differential spectra to three different 
models introduced above: the OPER model, the modified Va¬ 
lencia and the Cao model and check whether a description of 
the data is possible without contribution from the dibaryon res¬ 
onance. The sensitivity of our data to the different mechanisms 
(double A(1232), A(1440) excitation, as well as higher lying 
resonances) can be studied. The total cross-section can be used 
for consistency checks of the different analysis, in comparison 
to the already measured pn cross-section, as dis¬ 


cussed in 11521 l53n . 

Our paper is organized as follows: The experimental proce¬ 
dure is described in Sec. |2l The features of the models used for 
the description of the data are introduced in Sec.|3] We present 
and discuss the experimental results and the comparison with 
the models in Sec.|4]and draw conclusions in Sec.|5] 


2. Experimental procedure 

The experimental data have been obtained using HADES 
ifs^ located at the GSI Helmholtzzentrum fiir Schwerionen- 
forschung in Darmstadt, Germany. HADES is a modern multi¬ 
purpose detector currently operating in the region of kinetic 
beam energies of up to 2 A GeV for nucleus-nucleus collisions. 
The detailed description of the set-up can be found in if^ . here 
we briefly summarize the main features relevant for the present 
analysis. 

HADES is divided into 6 identical sectors defined by the su¬ 
perconducting coils producing the toroidal geometry magnetic 
field. The spectrometer has 85% of azimuthal acceptance and 
covers polar angles from 18° to 85° measured relatively the 
beam direction. Each sector of the spectrometer contains a Ring 
Imaging Cherenkov detector (RICH) operating in the magnetic 
field-free region, 4 planes of the Multi-wire Drift Chambers 
(MDC) located before and after the magnetic held region, two 
plastic scintillator walls for the polar angles larger (TOE) and 
smaller (TOEINO) than 45°, respectively, and an electromag¬ 
netic cascade detector (Pre-Shower) behind TOEINO for par¬ 
ticle identihcation. A two-stage trigger system is employed to 
select events within a predehned charged particle multiplicity 
interval, as well as the electron candidates. The investigation of 
the quasi-free np reactions with the deuteron beam is performed 


by using a Eorward Wall (EW) scintillator hodoscope by regis¬ 
tering the spectator protons. The EW is an array which consists 
of nearly 300 scintillating cells each 2.54 cm thick. During the 
dp experiment the EW was located 7 meters downstream the 
target covering polar angles from 0.33° up to 7.17°. A Monte 
Carlo simulation for deuteron-proton breakup has shown that 
approximately 90% of all spectator protons are within the EW 
acceptance 1171]. 

In the present experiment, the deuteron beam with intensity 
up to 10^ particles/s and 1.25 GeV/u of kinetic energy was di¬ 
rected onto a 5 cm long liquid-hydrogen target of 1% interac¬ 
tion probability. The analysis of the deuteron induced quasi- 
free np reactions was based on the first-level triggered events, 
with the required charged particle multiplicity of at least three 
in TOE and TOEINO and a signal in the EW hodoscope. The 
information from RICH, TOE/TOEINO and Pre-Shower detec¬ 
tors was used to select and discriminate the hadrons against the 
leptons. The momentum reconstruction was carried out by mea¬ 
suring the deflection angle of the particle trajectories derived 
from the four hit positions in the MDC planes (two before and 
two after the magnetic field zone) using the Runge-Kutta algo¬ 
rithm i^ . The achieved momentum resolution was 2-3% for 
protons and pions depending on their momentum and scattering 
angle. 

Quasi-free np interactions have been selected by the detec¬ 
tion of the proton spectators with scattering angles <2° and mo¬ 
menta between 1.7 GeV/c and 2.3 GeV/c reconstructed from 
the time-of-flight measurement in EW. The latter window is 
centered on half of the beam momentum and has a half width 
equals to three times the width of the momentum resolution. 
The np npn^n^ channel selection was based on identifica¬ 
tion of three detected hadrons where one of them has the nega¬ 
tive momentum polarity. The particle identification was based 
on event hypotheses where any of the three selected hadrons has 
been used as the reference particle. The reference particle time- 
of-flight was calculated using the reconstructed momentum and 
trajectory length. The velocities of the other two particles were 
then deduced, using only the time-of-flight difference with the 
reference particle. The applicability of this time-of-flight re¬ 
construction algorithm was checked in a dedicated experiment 
with a low beam intensity by means of the START detector as 
discussed in The correlations between the velocity and 

reconstructed momentum for all three particles were taken into 
consideration to reject the wrong hypotheses. The correlation 
between the energy losses in MDGs and momentum was ad¬ 
ditionally applied for the final selection of the protons and 
Einally, the total number of selected events corresponding to the 
np —> npji^n^ channel was ~ 8 ■ 10^. 


The normalization of the experimental yield has been per¬ 
formed using the simultaneously measured quasi-elastic pp 
scattering yield 05411 . The selection of pp elastic events was 
based on relations between the polar angles 6i and O 2 and az¬ 
imuthal angles 4>\ and (p 2 of both protons due to the momentum 
conservation in elastic scattering : 


= 180°, 


( 1 ) 
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been determined by means of full GEANT simulations, with 
pions and protons processed through the detector and recon¬ 
structed by tracking, particle identification and selection algo¬ 
rithms with the same package as it was done for real events. The 
acceptance of HADES for the np ^ npji^ji" reaction is about 
6% in comparison with the full phase-space. The resulting de¬ 
tection and reconstruction efficiency is typically about 90% for 
protons and pions. These efficiencies have a very smooth be¬ 
havior as a function of the different observables shown in the 
following. Therefore, we only consider a global systematic un¬ 
certainty which will affect the normalization of our data, but not 
the shape of the spectra. This effect is partially taken into ac¬ 
count by the normalization to the elastic data, and we estimate 
the residual uncertainty to 2%. 


Figure 1: Squared missing mass distributions of the pn^n^ system for the np —> 
npn*n^ reaction at 1.25 GeV. The experimental data are shown by the solid 
symbols. The shaded area displays the phase-space distributions at 1.25 GeV, 
corrected for the energy dependence of total cross section taken from l20h . 
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where jcm is the Lorentz factor of the center-of-mass (CM) 
system. The elastic events were approximately selected by an 
elliptic cut in the \(j)\ - (p 2 \ vs. tan0i ■ tan 02 ) plane with semi 
axes corresponding to 3cr cut for each variable. The quasi-free 
pp elastic data measured by HADES in the angular range 46° - 
134° in the CM system have been corrected for the efficiency 
and acceptance by using a Pluto ll55ll56l] simulation which uses 
the data from 15411 as an input and takes into account the proton 
momentum distribution in the deuteron and the dependence of 
the cross section on the pp center of mass energy dependence 
of the cross section. The resulting factor taken to normalize the 
differential cross section has a precision of about 4% at 1.25 
GeV. The details of the normalization procedure at 1.25 GeV 
are discussed in 11571] . 

The purity of the np —» npn^n^ channel selection is demon¬ 
strated by the squared missing mass distribution of pn^n^ 
shown as solid circles in Eig. [1] which peaks close to the neu¬ 
tron mass. The shaded area displays the result of a simulation 
of the quasi-free np —» npji^n^ reaction, where the neutron mo¬ 
mentum distribution in the deuteron is taken into account using 
the Paris potential ll^ and the production in the n + p 
reaction is treated using phase-distributions and considering a 
rise of the cross section with the np center-of-mass energy ac¬ 
cording to i^ . 

The direct comparison of theoretical and efficiency cor¬ 
rected experimental distributions has been performed inside the 
HADES acceptance using dedicated filters, as will be shown in 
the following. Eor this purpose, the acceptances and efficiencies 
for different particles (i.e. pions and protons) were separately 
tabulated in matrices as functions of the momentum, azimuthal 
and polar angles. The acceptance matrices describe only the 
HADES fiducial volume and can be applied as a filter to the 
events generated by the models. The matrix coefficients have 


3. Model features 

Eor the interpretation and discussion of our results in sec. |4] 
below, we used the modified Valencia 11331] . the Cao ll46ll and 
OPER models ]l43ll . The Cao and modified Valencia models are 
effective lagrangian models taking into account different res¬ 
onant and non-resonant graphs. The original Valencia model 
includes only V(1440) and A(1232) resonances while the 
Cao model considers all known baryonic resonances with mass 
up to 1.72 GeV, neglecting the interferences between the differ¬ 
ent contributions, included in the Valencia model. The practical 
differences in our energy range are lower contributions for both 
A and V(1440) by about 30% in the Valencia model, which are 
partially compensated by constructive interferences. In com¬ 
parison with the original Valencia model, four main changes 
have been applied in the ’’modified” Valencia model. Eirst, the 
ratio R= r(V(1440) — > A7r)/r(V(1440) Ncr) of the branch¬ 
ing ratios of the Roper resonance towards 27rN via An or Ncr has 
been reduced from 4 in the initial model to 1 . The initial value 
corresponded to PDG ifs^ estimates prior to 2012, and a new 
value was deduced from an analysis of the opening angle 
and invariant mass distributions obtained in the pp ppn^^n^^ 
reaction below 900 MeV assuming the dominance of the 
Roper excitation. This new value is also in agreement with a 
recent Partial Wave Analysis M]. In the meantime, the 
PDG limits for the ratio have also been changed and are now 
between 1 and 3. This change seems therefore fully consistent. 
The Cao model uses a value 2 for this ratio, hence favoring 
V(1440) —> Att by a factor 2 with respect to V(1440) —» Ncr 
decay, while they have equal weights in the modified Valencia 
model. 

The second change is a readjustment of the V(1440) strength 
in the modified Valencia model according to the isospin decom¬ 
position of the two-pion production channels in the pp reaction 
at different energies between 0.775 GeV and 1.36 GeV [31]. 
The maximum reduction is obtained exactly for our incident 
energy of 1.25 GeV and amounts to a factor 2. This brings the 
V(1440) contribution to be much smaller than the double A ex¬ 
citation in the modified Valencia model and makes a significant 
difference with the Cao model, where the V(1440) contribution 
corresponds to about 40 % of the total pn —» pnn^n" cross 
section at an incident energy of 1.25 GeV. 
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The third modification is driven by the shape of the in¬ 
variant mass distribution measured in pp —> ppji^Ti^ at an inci¬ 
dent energy larger than 1 GeV, where the double A mechanism 
dominates 1^ . It consists in a reduction by a factor 12, and 
a change of the sign of the p exchange contribution in the A 
excitation mechanism. The latter indeed induced a two-hump 
structure in the 7r‘’;7r*' invariant mass distribution, which was not 
seen, neither in the data nor in the Cao model which has a 
very small p exchange contribution. 

The fourth modification consists in the introduction of the 
A(1600) ^ A(1232)7r contribution to improve the description 
of the pp —» nmr^TT^ reaction. In this channel, the N* decay 
does not contribute, and the cross section is larger by more than 
a factor 2 than expected for the double A(1232) contribution 
only. The WASA collaboration explained this feature by a large 
contribution of A(1600) excitation, with a destructive interfer¬ 
ence with the double A(1232) contribution. After introducing 
this new contribution, the pp ^ nnn^n'^ cross section could be 
reproduced and the description of the differential spectra, espe¬ 
cially the opening angle and invariant mass distributions 
measured for this channel at 1.1 GeV improved signifi¬ 
cantly. The A(1600) —» A(1232)7r process was implemented in 
a similar way as A(1232) ^ A(1232);7r, which was already taken 
into account in the Valencia model. However, such a large con¬ 
tribution of the A(1600) resonance at such low energies is sur¬ 
prising. In the Cao model, it is taken into account, but starts to 
contribute significantly to the pp ^ nnn^TT^ reaction for ener¬ 
gies larger than 1.6 GeV. It has to be noted that, while the orig¬ 
inal Valencia model underestimates the pp cross 

section by more than a factor 2, the Cao model provides good 
pred ictions for both total cross sections and differential spectra 
IIThll . This is related to the larger double A contribution in this 
model. However, such a large double A contribution leads to a 
significant overestimate of the pp —» ppn'^n^ cross section. In 
both models, the nucleon pole terms are found to be very small, 
although they have a much larger relative contribution in the 
pp —> nnn^n^ than in other channels. 

As mentioned in the introduction, the OPER model is based 
on a Reggeized tt -exchange model and uses on-shell ampli¬ 
tudes of the elastic nN scattering and of the inelastic nN 
nnN reaction, with form factors and propagators taking into 
account the off-shellness of the exchanged pion. The elastic 
amplitudes are taken from a Partial Wave Analysis and the in¬ 
elastic ones are deduced from a parametrization obtained in 
the framework of the Generalized Isobar Model ^]. This 
model provided a good description of differential spectra mea¬ 
sured at Dubna in the np npn^n" reaction at 5.2 GeV/c and 
in pp —> ppn'^n^ at 7.23 GeV/c ll^ . The One-Baryon Ex¬ 
change (OBE) diagrams were introduced, as described in 13811 . 
to improve the description of the np —> npn^n^ spectra below 
3 GeV/c. Such diagrams correspond in Lagrangian models to 
the ’’pre-emission contribution”. Very recently, the contribu¬ 
tion of the so called ’’hanged” diagrams, due to the two-pion 
production from the exchange pion line, was also considered. 
Eor this, amplitudes for nn scattering were used. This addition 
was shown to improve the description of the nn invariant mass 
spectra in the low mass region 11431] . Such graphs are taken into 


account in the Valencia model, but are neglected in the Cao 
model. Table [T]displays the main resonant contributions in the 
three models at an incident neutron energy of 1.25 GeV, disre¬ 
garding the contribution via interference terms. This makes a 
big difference with the OPER model, where the OBE diagrams, 
where at least one pion is not emitted by a resonance, amount to 
41% of the total yield. Although these numbers do not take into 
account interference effects, they already point to major differ¬ 
ences between the models which will be investigated further 
when comparing to our data. 



Cao 

mod. Valencia 

OPER 

A(1232)A(1232) 

47.0 % 

60% 

38.0 % 

At(1440) ^ A(1232)7r 

23.0 % 

2.1 % 

4.5 % 

At(1440) ^ Ncr 

20.0 % 

8 .2% 

0 .2% 

A(1600)^ A(1232)7r 

3.0% 

21 .0% 

4.5 % 


Table 1: Main contributions in the Cao (H, ’’modified Valencia” (bSi and 
operB models for the np —> npn^rr for a neutron incident energy of 1.25 
GeV. 


To take into account the momentum distribution of the neu¬ 
tron inside the deuteron, the Paris if^ deuteron wave functions 
has been used for the phase-space calculation, the Hulthen ll63ll 
wave function for modified Valencia ll^ and CD-Bonn 1 64 ] for 
OPER 1^] and Cao models. No significant difference is 
expected from these different inputs. As mentioned in the pre¬ 
vious section, events have been generated according to distribu¬ 
tions from these models and filtered by the HADES acceptance. 


4. Experimental results and comparison with models 


We first concentrate on the comparison of the shapes of the 
theoretical and experimental distributions and will discuss the 
cross sections at the end of this section. The invariant mass 
and angular distributions for the «/? —> npji^n^ reaction at 1.25 
GeV measured inside the HADES acceptance and corrected for 
the reconstruction efficiency are presented by solid circles in 
Eigs.|2]l3]and|4]l5]respectively. Only statistical errors are shown. 
As explained above, the global uncertainty taking into account 
both normalization and efficiency correction is of the order of 
5%, with no significant dependence on the different observ¬ 
ables. 

The predictions from the modified Valencia |33], Cao iS 
and OPER models ll^ have been normalized therefore to the 
total experimental yield and are presented as long-dashed, 
dashed and solid lines, respectively, in Eigs. |2] and |4] The 
hatched areas correspond to phase-space distributions obtained 
in the same way. In addition, the different contributions in the 
OPER and modified Valencia models are shown in Eigs.[3]and|5] 
in comparison to selected experimental distributions exhibited 
in Eigs. |2] and 0] 


4.1. Invariant mass distributions 

The experimental pjt^, pn^ and pn^n~ invariant mass distri¬ 
butions (panels b), c) and d) in Eig. |2]i are all shifted to higher 
masses in comparison with the phase-space calculations. In 
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Figure 2: Distributions of the (a), p-K~ (b), (c) and pn^TC (d) in¬ 

variant masses for the np npn^n~ reaction at 1.25 GeV. The experimental 
data are shown by solid symbols. The theoretical predictions within HADES 
acceptance from OPER f4^ . Cao da and modified Valencia models are 
given by the solid, dashed and long-dashed curves, respectively. The shaded 
areas show the phase-space distributions. 


particular, the invariant mass distribution of the subsys¬ 
tem (panel c) in Fig. |2]i shows a pronounced resonant behavior 
with a position of the maximum in the experimental distribu¬ 
tion roughly corresponding to the mass. This distribution 
is well described by the different models and deviates signih- 
cantly from the phase space distribution. Such a resonant be¬ 
havior of the pTT^ invariant mass distribution is expected for the 
double A excitation process, since the isospin factors favor the 
excitation of the A^^A“ contribution with respect to the A^A® 
by a weight of 8/5. A resonance in the pn^ system is also ex¬ 
pected for the A(1440)^ decay into but with a ’’smear¬ 

ing” due to the excitation of the A(1440)° and subsequent decay 
into A TT^ which has the same probability. In addition, the dom¬ 
inance with respect to the other isospin channels is also lower 
than in the case of the double A excitation. The situation is sim¬ 
ilar for the A(1600) case, with an even lower relative weight for 
the decay involving the A(1232)^^ excitation. These statements 
are confirmed by the behavior of the A(1440) and A(1600) con¬ 
tributions which are displayed in Fig. |2l below in the case of 
the modified Valencia model. The overall contribution of the 
double A excitation is larger in the modified Valencia model 
than in the Cao model. This is probably the reason why the lat¬ 
ter presents a slightly broader pn^ invariant mass distribution. 
Surprisingly, the pn^ invariant mass of the double A(1232) con¬ 
tribution in the OPER model is slightly broader than in the Va¬ 
lencia model, as can be seen in Fig.|3};. The global distribution 
which is further broadened by the OBE contributions, is similar 
to the one in the Cao model with an additional shift to higher 
energies (see Eig.|2]r). 

The deviations from phase-space distributions are less spec- 



Figure 3: Experimental data (full dots) for np —* npn*!T^ at 1.25 GeV are com¬ 
pared to the total yield (solid curves) for the OPER model l43ll (left column) 
and modified Valencia model |3^ (right column). In each case, the AA 
(long-dashed), V(1400) (long dash-dotted) and A(1600) (short-dashed) contri¬ 
butions are shown, a) and b): invariant mass , c) and d): pn^ invariant 

mass, e) and f): pn^n^ invariant mass distributions. 


tacular for the pn^ and pn^n" invariant mass distributions (pan¬ 
els b) and d) in Eig. |2]l. No clear resonance behavior is in¬ 
deed expected for the pji^ system, since the A(1232)° is disfa¬ 
vored by isospin with respect to A^^ in all channels (double A, 
N ^A(1232);r, and the A(1600)^A(1232)7r). This also holds 
for the pn^n^ system, due to the strong double A contribution. 
In addition, the A(1440)° excitation, which is as probable as 
the A(1440)^, is favored by acceptance. This can be checked 
in Figs. [3}; and[3}l below, where the A(1440) contributions are 
indeed shifted to low pji^n" invariant mass. For the A(1600) 
excitation, the positive charge state, decaying into pn^n^, is 
favored by acceptance. This is probably the reason why the 
pn*n~ invariant mass distribution for this contribution is shifted 
in the models to the high-mass part of the available phase space. 
The OPER and Cao models give similar predictions for the pn^ 
and p 7 T^ 7 T'~ observables and achieve a reasonable description of 
the pn~ but predict a too narrow pn^7t~ invariant mass distribu¬ 
tion. In the case of the Valencia model, both the pn^ and pn^n^ 
are overestimated on the high energy side. This points to a too 
large A(1600) contribution, as will be discussed below in some 
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details. 

The three models present the largest deviations for the 
invariant mass distribution. The sensitivity of this distribu¬ 
tion to the two-pion production mechanism in NN —> NNnn 
reactions has already been demonstrated in previous work 
1 3 S 12 13l- In fact, the puzzling enhancement observed 
in the low-energy part of the invariant mass spectra of two pi- 
ons produced in pn or pd fusion reaction (the so-called ABC 
effect) was also a manifestation of this sensitivity and triggered 
the onset of double A excitation models llhsll . The OPER model 
is not too far from the data, but has a slightly too flat shape, 
while the Cao and Valencia models fail strongly to describe the 
shape of this distribution. The double-hump structure which 
can be seen in the predictions of the Cao model is completely 
absent in our data. A very similar trend was already observed 
for the comparison of the Cao model predictions to the spectra 
measured in the jw ppn^n^ and pp ppn^n^ reactions 

above 1.1 GeV 11461] . Since this double-hump structure is due to 
the V(1440)^ Att process, these results indicate that this con¬ 
tribution is too large in the Cao model. The OPER model gives 
the best description of the low mass part of this distribution, 
which seems to be mainly due to the addition of the OBE and 
’’hanged” diagrams. However, the double A contribution is also 
much broader than in the case of the modified Valencia model 
(see Eigs. and b below). In the modified Valencia model, 
the excess in the region above 350 MeV/c^ is due to both the 
A(1600) and the Roper contributions and their interference with 
the double A(1232) contribution. In particular, the interference 
between the Roper and double A(1232) contribution is small in 
overall, but has a significant constructive effect of the level of 
10-15% in this region. Our data would be better described by a 
change of sign of this interference together with a reduction of 
the A(1600) contribution. 



Figure 4: Same as Fig. |2] but for angular distributions: a) - opening angle of 
in the np rest frame, b) - polar angle of in the np rest frame, c) 

- polar angle of pn^ in the np rest frame, d) - polar angle of pn* in the np 
rest frame, e) - polar angle of n^ in the pn^ Gottfried-Jackson frame, f) - polar 
angle of n* in the pn* Gottfried-Jackson frame. 


4.2. Angular distributions 

Eigure |4] exhibits the angular distributions for the np —> 
npn^n^ reaction at 1.25 GeV. Panels a), b), c), d) and e), /) cor¬ 
respond to the distributions of the opening angle of polar 
angles of pn^, pn^, pn^ji^ subsystems in the np CM system 
and polar angles of and ji^ in the pn^ and pn^ Gottfried- 
Jackson frames, respectively. The CM frame was defined as¬ 
suming the neutron at rest in the deuteron and the angles in the 
Gottfried-Jackson frame are defined with respect to the beam 
direction. 

The excess of the measured yields at small opening an¬ 
gles (cos ~ 1) with respect to the predictions of all models 
observed in Pig. |4] reflects the enhancement at small in¬ 
variant mass which was observed in Pig. |2^. Both variables 
are indeed strongly correlated. The large asymmetry of the 
distribution in the case of the Cao model and its steep peak¬ 
ing for back to back ji^n^ emission (cosd™- = -1) are re¬ 
lated to the high-mass structure in Pig. |2^, which, as already 
mentioned, is due to the V(1440) —>A(1232) ji decay. As ex¬ 
plained in 1 ^ . this effect is produced by the double p-wave 
decay which gives an amplitude with a cos^ behavior. Por 
the other models, the backward peaking is consistent with the 


data. It can be checked in Pigs. |5^ and b that the distribu¬ 
tions are rather smooth for all contributions. In the case of 
the modified Valencia model, the roughly symmetric distribu¬ 
tion results from the opposite trends of the double A(1232) and 
A(1600) which have yields respectively forward and backward 
peaked in the HADES acceptance. Por the OPER model, the 
anisotropy is mainly due to the OBE terms, all other contri¬ 
butions being very flat. Since the V(1440) contribution in the 
OPER model is dominated by the Att decay, we were expecting 
a larger anisotropy, following the arguments given above. In the 
case of the modified Valencia model, the V(1440) A(1232)7r 

decay is suppressed with respect to the crN decay (see TablellJ. 
It is therefore not surprising that the V(1440) contribution does 
not present a strong anisotropy. The difference of the shapes of 
the cos dS*!- distributions obtained for the double A contribu- 
tions in the modified Valencia and OPER models is consistent 
with the already mentioned different behavior of the in¬ 
variant masses (Pigs. |5^ and b). 

The and 0^^ angular distributions shown in pan¬ 

els b) and c) and d) panels in Pig. |4] respectively, present a sig¬ 
nificant forward/backward asymmetry. This asymmetry reflects 
the fact that protons are emitted preferentially backward in the 
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Figure 5: Same as Fig. [^for opening angle of ti^-k in the np rest frame ( a) 
and b)) and polar angle of pn* system in the np rest frame ( c) and d)). 


CM, which is mainly due to the strong contribution. 

This process is more likely for small four-momentum transfers 
between the target proton and the A^^ and the beam neutron 
and the A“ respectively, hence the peaking of the protons at 
backward angles. Such a strong asymmetry is not present for 
the other contributions. However, the effect is also enhanced by 
the acceptance, as demonstrated by the distribution of the pure 
phase space. Due to the lacking detection capabilities for labo¬ 
ratory polar angles lower than 18°, the detection of the proton 
is indeed more likely at backward center-of-mass angles. The 
Cao model gives the best overall description of the three dis¬ 
tributions. The modified Valencia model is also doing rather 
well, however forward/backward asymmetry is smaller than in 
experimental data. A much worse description of the slopes of 
(f^ and (f^ is obtained in OPER model. As shown 
in Figs. 1^ and |2l, the double A contribution is much more 
strongly backward peaked than in the modified Valencia model. 
However, the OBE and ’’hanged” contributions, which are not 
displayed in the picture, are also responsible for the too steep 
total angular distribution. One of the possible reasons of such 
a deviation can be the omission of the interference between the 
’’hanged” and other diagrams ll43l] . The striking difference with 
the double A contribution in the modified Valencia model is 
probably due to the much lower cut-off parameters in the ver¬ 
tex form factors which induce a much steeper four-momentum 
transfer dependence. Worth to note is the strong difference be¬ 
tween the and 6*'^ distributions, in contrast with the phase 
space distribution, which allows to appreciate the small effect 
of the different and acceptances. Again, the modified 
Valencia model gives the best simultaneous estimate of these 
distributions, which we take as a hint that this model correctly 
describes the different isospin configurations of the two-pion 
production mechanisms. It can also be noted that, similarly to 


the TT^TT invariant mass distributions, the description of the 0^“ 
distribution could be improved by a reduction of the A(1600) 
contribution. The distributions of the pion angles (GJ) and 
(GJ) are shown in panels e) and /) in Fig. |4] respectively. 
The quantity 0y denotes the angle between the thee-momentum 
of particle j and the direction of the beam in the center-of-mass 
of particles i and j, taken relatively to the direction of the beam 
particle in this reference frame. These distributions also present 
striking differences for and ttA It can be noted that the 
0^^ (GJ) distribution is very well reproduced, especially by the 
OPER and modified Valencia model, while none of the models 
predict the observed enhancement for backward 0^f (GJ). 

To summarize the analysis of these distributions, the Valencia 
model provides a much better description than the Cao model. 
Our analysis therefore validates the changes introduced by the 
WAS A collaboration in the original Valencia model, except per¬ 
haps for the A(1600) contribution, which seems to be too large. 
The OPER model, which is based on a very different approach, 
gives also a good description of the data. In particular, it fits bet¬ 
ter the 7r^7r“ invariant masses than the modified Valencia model. 
However, the predictions are worse for the CM angular distri¬ 
butions of the pn^, pn wA pn^n" subsystems. 



Figure 6: FIADES measurement for the quasi-free np —> npn^n^ reaction us¬ 
ing a deuterium beam at 1.25 GeV/nucleon (full dot) compared to world data 
shown by various symbols. The horizontal error bars indicate the spread of 
the neutron momentum in the ditferent measurements. The full and short dash- 
dotted curves display respectively the "Bistricky parametrization” used for the 
OPER model 1^ normalization and the predictions of the modified Valencia 
model Is^l . The long dash-dotted curve is the estimate from IsSl for the contri¬ 
bution of the dibaryon resonance. The dashed curve is the sum of the modified 
Valencia model and dibaryon resonance contributions. 


4.3. Absolute cross sections 

We come now to the discussions of the absolute yields. The 
measured differential cross section integrated over the HADES 
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acceptance is 34.9 +1.5 pb. The model 114611 and Il43h predic¬ 
tions, respectively 72.4 and 86 /ib are larger by a factor more 
than 2, while the modified Valencia model 1331] is in better 
agreement, with a cross section of 26.4 pb, i.e., about 30% 
lower than the experimental value. Acceptance corrections 
were calculated using the modified Valencia model 13311 and 
OPER model which give a reasonable description of the dif¬ 
ferential distributions, as shown in Figs. |2]|5] In practice, eight 
different differential spectra, corresponding to the distributions 
in Mjr+rr, M,„r, M cos™ ,cos™ cos0™ and 

n n 1 pn 1 pn i pn n i 6 „+„- ’ pn^n ’ pr 

cos 6^^ (Figs. |2] and |4]a-d) were corrected for acceptance us¬ 
ing both models. The differences between the integrated yields 
obtained for these different distributions were used to deter¬ 
mine the systematic errors. Following this procedure, we ob¬ 
tain for the total cross section 0.65+0.03 mb using the modified 
Valencia model and 0.795+0.040 mb using the OPER model. 
Our final estimate, taking into account these two values is cr- 
0.722+0.108 mb. This cross section is averaged over neutron 
energies accessible in the quasi-free np —> npn^n^ reaction 
at a deuteron beam energy of 1.25 GeV/nucleon. According 
to the fit of the missing mass (Fig. [T]), taking into account the 
neutron energy and momentum distribution in the deuteron, the 
average neutron energy for this measurement is 1.273+ 0.063 
GeV. Our data point is shown together with the world data in 
Fig. |6l In contrast to previous plots which can be found in the 
literature, we take into account the spread of the neutron mo¬ 
mentum in the different measurements. The Cao model over¬ 
estimates our point by a factor 2.4, with a prediction of 1.730 
mb, which confirms that this model does not reproduce satis¬ 
factorily the np —> npn^n" reaction in our energy range. The 
OPER model does not provide cross sections and was normal¬ 
ized to the ’’Bistricky parametrization” (shown as a black 
curve in Fig.|6]), which resulted from a simple interpolation be¬ 
tween measurements over a wide energy range up to 2.2 GeV. 
This parametrization largely overestimates the NIMROD mea¬ 
surement lll4ll . which was obtained at an incident energy higher 
by about 120 MeV than our experiment. On the other hand, it 
provided a good prediction for the measurement at KEK in 
which was obtained in the meantime. The latter is of special in¬ 
terest for the present study, since it was obtained for an incident 
energy only 70 MeV lower than our experiment. Our extrapo¬ 
lated cross section is approximately a factor 2.6 lower than the 
’’Bistricky parametrization” that has a value of 1.88 mb at 1.25 
GeV. This result is in contradiction with a smoothly increasing 
cross section as a function of the incident energy. This, as we 
will see in the following, could be expected in the presence of 
a resonance at lower energies. Our measurement is however 
hardly compatible with the KEK data (cr = 1.25+ 0.05 mb at 
T„=1.17 GeV). A decrease of a factor 2 in such a small energy 
range is indeed difhcult to explain. On the other hand, the mod¬ 
ified Valencia gives a prediction of 590 ph, i.e. much closer 
to our value. This model also reasonably reproduces not only 
the differential distributions for the np ^ npji^n^ channel, as 
shown in this paper, but also the total cross sections and differ¬ 
ential distributions for pp ppn^n^ 1321] and pp nnn^n^ 
channels measured by WAS A below 1.4 GeV. However, 


this model does not describe satisfactorily the general excita¬ 
tion function of the np —> npn^n" reaction. The underestima¬ 
tion of the data at higher energies might well be due to the lack 
of higher lying resonances, such as V(1520) and V(1535) in 
the model. However, another explanation has to be found for 
the underestimation at lower energies. 

In the hypothesis of a dibaryon resonance with a mass around 
2.38 GeV, as claimed by the WASA collaboration, a structure is 
expected at a neutron energy around 1.13 GeV. Using the cross 
sections for pn d* extracted by the WASA collab¬ 

oration, Albaladejo et al. if^ have estimated the resonant cross 
section for the np —> np;7r^;7r“ channel, i.e., pn —> —» pnn^n~ 
(see dash-dotted curve in Fig. |6ll. They concluded on the dif¬ 
ficulty to reconcile the existing data with the resonance hy¬ 
pothesis, taking into account the large non-resonant contribu¬ 
tion which has necessarily to be added. However, this con¬ 
clusion is, to our opinion, strongly biased by the KEK point 
which is already beyond the resonance and favors a large non¬ 
resonant contribution. In this respect, our measurement, which 
provides a much smaller contribution, is in better agreement 
with the dibaryon resonance hypothesis. The expected reso¬ 
nant contribution at our energy is about 0.2 mb, less than one 
third of our measured cross section, and is compatible with a 
non-resonant contribution of the order of 0.5 mb as predicted 
by the modified Valencia model. According to this model, the 
non-resonant contribution at the peak of the resonance is only 
of the order of 0.2 mb. The existing cross section measure¬ 
ments in this region are therefore not inconsistent with a reso¬ 
nant cross-section of about 0.7 mb, as deduced from the WASA 
result in the pn dn^n^ channel. To illustrate that, we show in 
Fig. 12 the excitation function obtained by adding the dibaryon 
contribution from ifs^ to the modified Valencia model predic¬ 
tion. Except for the KEK point, the result describes within 
20 % the cross sections obtained for incident energies between 1 
and 1.3 GeV, which is the energy range of the pp—>pp7r‘'.7r° and 
pp—>nn;7r+7r^ data used to adjust the modified Valencia model. 
In particular, it is in perfect agreement with our measurement. 
The poorer description outside this energy range could proba¬ 
bly be reduced by further adjustments of the Valencia model. 
No conclusion can however be drawn from this very crude cal¬ 
culation, due to the unknown effect of interferences between 
t-channel and s-channel processes. In addition, it would be nec¬ 
essary to include the decay of the resonance into a pn pair with 
other quantum numbers than the deuteron. According to the 
analysis in this contribution is estimated to be of the order 
of 0.1 mb. The effect of the resonance contribution on the differ¬ 
ential distributions and in particular to the pji^n" invariant mass 
is also an open question. This can only be made using a full 
model, including in a consistent way the t-channel processes, 
based on the ’’modified Valencia” model and the resonant s- 
channel although a phenomenological approach was presented 
for the analysis of np ^ npifiifi in S. 

4.4. Acceptance corrected distributions 

In order to allow for a more direct comparison of our data to 
differential cross sections for quasi-free np —> npn^n" lf3^ and 
np —> npifiifi ifj^ reactions measured in other experiments and 
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Figure 7: The acceptance corrected distributions of the invariant mass 

(a) and opening angle of n'^n~ in the np rest frame (b) for the np np7T^n~ 
reaction at 1.25 GeV. The experimental data are shown by solid symbols. The 
theoretical predictions in 4;r region from OPER O . Cao (H and modified 
Valencia models are given by the solid, dashed and long-dashed curves, 
respectively. The shaded areas show the phase-space distributions. 


possibly facilitate the comparison to other potential models key 
distributions have been corrected for acceptance. The correc¬ 
tion has been done by using modified Valencia lO and OPER 
ll43l] models. Fig. [T] presents the extrapolated distributions of 
;7r^7r“ invariant mass and opening angle between pions in CM 
frame for np —> npn^n^ reaction at 1.25 GeV. The errors on 
extrapolated data take into account the difference between ac¬ 
ceptance correction coefficients for modified Valencia 03311 and 
OPER lldlll models. The predictions of the different models 
in 47r are also displayed in Fig. Q using the same normalisa¬ 
tion as for the comparison inside HADES acceptance (Fig. 

- |5]l. The different yields therefore illustrate the different ac¬ 
ceptance correction factors. The same features as described in 
Sec. 14. II and Sec. 14.21 for the results inside acceptance can be 


observed for the results extrapolated in 47r. Indeed the OPER 
and modified Valencia ll^ models give the best overall 
description of the shape of both distributions. The OPER ll43l] 
and in a lesser extent the modihed Valencia models give 
also good results for the invariant mass and opening an¬ 
gle distributions measured in the np —» npji^n^ at 1. and 1.5 
GeV. The results obtained at the three energies allow to follow 
the evolution of the mechanisms, characterized in particular by 
the important role of non-resonant contributions at 1 GeV and 
the dominance of double Delta excitation at 1.5 GeV. As for the 
np npn^^TT^^ reaction at 1.13 GeV iIt^ . where a dominance 
of the d* contribution is observed, the shape of the nn invariant 
mass is similar to our measurement, which confirms the rather 
low sensitivity of this observable to the s-channel resonant con¬ 
tribution, as stressed by the authors. 


5. Conclusion 

We have presented an analysis of high-statistics differential 
distributions obtained in the np ^ npn^n^ reaction, using a 
deuterium beam at 1.25 GeV/nucleon and the HADES experi¬ 
mental set-up at GSI. In this channel, in addition to t-channel 
processes, mainly due to the Roper and double A(1232) excita¬ 
tions, the contribution of a s-channel process, with an interme¬ 
diate dibaryon, as observed in the reaction pn dn*n'~ might 
be present. In this paper, we focused on the compatibility of our 


data with models including t-channel processes only. The ex¬ 
perimental distributions have therefore been compared with the 
predictions of two Lagrangian models, the Cao model 11461] and 
the ’’modified Valencia” model lf33l] and of the OPER model 
iilll . based on one-pion and one-baryon exchanges. This anal¬ 
ysis demonstrates the high sensitivity of the differential distri¬ 
butions to the different components of the two-pion production 
mechanism. The modified Valencia model gives a much bet¬ 
ter description of the shapes of the distributions than the Cao 
model. This good result supports the changes introduced by the 
WASA collaboration to describe the total and differential cross 
sections in pp —> ppn^n^ and pp —> nnTr^rr^ in the same energy 
range, especially the reduction of the V(1440) —>A(1232) n 
contribution. The double A contribution seems to be rather well 
described by the model, which is important since this process is 
expected to contribute to dilepton production in the pn reaction 
ifsl] . However, some discrepancies are observed in our channel, 
like a missing strength at small invariant masses and for¬ 
ward angles in the model. These discrepancies could be 
probably reduced by further adjusting the A(1600) and Roper 
contributions. However, such changes should also be validated 
by an analysis of two-pion production in other isospin channels. 
The OPER model gives also rather good results, except for the 
polar angle of the pn"^ system in the np rest frame, with a lower 
relative contribution of the double A(1232) contribution, and a 
significant contribution of one-baryon exchange graphs. Both 
models have been used to calculate acceptance corrections, and 
obtain an estimate of the total cross section of the np ^ npn^n^ 
reaction at 1.25 GeV beam energy. This measurement is im¬ 
portant due to the scarce existing measurement in the relevant 
energy range. The found rather low value of the cross section 
is not inconsistent with a resonant structure at low energies, as 
expected in presence of the dibaryon resonance, with mass M~ 
2.38 GeV and width T ~ 70 MeV reported by the WASA collab¬ 
oration. However, our measurement is hardly compatible with 
the measurement performed at KEK at an incident neutron en¬ 
ergy lower by 70 MeV. The modified Valencia model, which has 
been now validated for different channels, underestimates the 
total cross section in our measurement by only 30%. Under the 
resonance peak, the underestimation of the data is much larger 
and is also compatible with the resonant hypothesis. However, 
the present situation, both from experimental and theoretical as¬ 
pects, is not clear enough to draw conclusions on the existence 
of the dibaryon resonance. On the experimental side, useful 
constraints in the present experiment could be obtained from 
the extraction of cross section at different neutron energies, the 
on-going analysis of the pp ppn^n^ and pn —> dn^n" chan¬ 
nels with HADES will allow for specific tests of the modihed 
Valencia model and the dibaryon resonance contribution, re¬ 
spectively. The experimental situation should therefore further 
become clearer in a near future. On the theory side, we would 
like to call for the development of a full model, including in 
a consistent way the t-channel processes, based on the modi¬ 
hed Valencia model and the s-channel processes including the 
dibaryon with above quoted parameters, which could provide a 
solid framework for the interpretation of the two-pion produc¬ 
tion data. 
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